Metaplectic Anyons, Majorana Zero Modes, and their Computational Power 



(N 

o 

> 
O 

r- 

(N 



Matthew B. Hastings, 12 Chetan Nayak, 2 - 3 and Zhenghan Wang 2 

'Duke University, Department of Physics, Durham, NC 27708, USA 
2 Station Q, Microsoft Research, Santa Barbara, CA 93106-6105 
3 Department of Physics, University of California, Santa Barbara, CA 93106 
(Dated: November 28, 2012) 

We introduce and study a class of any on models that are a natural generalization of Ising anyons and Majorana 
fermion zero modes. These models combine an Ising anyon sector with a sector associated with SO(m)2 
Chern-Simons theory. We show how they can arise in a simple scenario for electron fractionalization and give a 
complete account of their quasiparticles types, fusion rules, and braiding. We show that the image of the braid 
group is finite for a collection of 2n fundamental quasiparticles and is a proper subgroup of the metaplectic 
representation of Sp(2n — 2, F m ) x H(2n — 2, F m ), where Sp(2n — 2, F m ) is the symplectic group over the 
finite field F m and H(2n — 2, F m ) is the extra special group (also called the (2n — 1) -dimensional Heisenberg 
group) over F m . Moreover, the braiding of fundamental quasiparticles can be efficiently simulated classically. 
However, computing the result of braiding a certain type of composite quasiparticle is #P-hard, although it 
is not universal for quantum computation because it has a finite braid group image. This a rare example of a 
topological phase that is not universal for quantum computation through braiding but nevertheless has #P-hard 
link invariants. We argue that our models are closely related to recent analyses finding non-Abelian anyonic 
properties for defects in quantum Hall systems, generalizing Majorana zero modes in quasi- ID systems. 
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PACS numbers: 



I. INTRODUCTION 

Majorana zero modes can occur in a wide variety of phys- 
ical systems linked by the common thread of chiral p-wave 
superconductivity and its analogs^—. They exhibit many 
(and, in some cases, nearly all) of the properties of Ising 
anyons and, therefore, may prove useful for fault-tolerant 
topological quantum information processin g 11 ' 20 . However, it 
is possible to classically simulate the braiding of Ising anyons 
efficientl y 21 ' 22 . Therefore, they are useful for quantum com- 
putation only if braiding is supplemented by measurement 
at intermediate stages of computations and by a ir/8 phase 
gate, in which case they are capable of universal quantum 
computation 23 . While it is likely that the former can be per- 
formed accurately, the latter appears difficult, although there 
are various interesting concrete proposals 2 ^—. Moreover, a 
non-topological implementation of the it/ 8 phase gate re- 
quires error correction, which entails significant overhead 26 . 
Therefore, physical systems supporting anyons that are capa- 
ble of universal quantum computation with braiding alon o 21,22 
(best-case scenario) or braiding and measurement 2 -^ (next- 
best scenario) would be a very attractive platform for quantum 
computation. 

In this paper, we introduce a sequence of topological 
phases of electrons which generalize physical models of Ising 
anyons. Suppose that an electron fractionalizes into a spinless 
neutral fermion ip and a charged spinful boson Z. Further, 
suppose that the spinless neutral fermion forms a p + ip paired 
superfluid state. If the bosons form a trivial gapped state, then 
the system is in the Ising anyon state, as in Kitaev's honey- 
comb lattice modeli (If the bosons condense, then the sys- 
tem is in a superconducting state which is a quasi-topological 
phase with some of the properties of Ising anyons^ 29 .) If the 
bosons form a spin-polarized fractional quantum Hall state, 
then the system is in the Moore-Read state 1 , the anti-Pfaffian 



stated, or a Bonderson-SlingerlancU^ state descended from 
one of these. But suppose, instead, that the bosons form a 
more complex topological phase of their own, T . Then the 
system will support quasiparticles that are combinations of 
those of the Ising topological quantum field theory (TQFT) 
and those of T, subject to the condition that they braid trivially 
with electrons. In the phases analyzed in this paper, T is asso- 
ciated with SO(m)2 Chern-Simons theory, where m = 3, 5, 7 
with, we believe, a generalization to any odd prime m. The 
50(m)2 TQFTs have several very interesting properties. All 
of these theories have a quasiparticle that is a boson. We 
identify this boson with Z through a non-Abelian analog of 
flux-attachment 3 ^—. In addition, these theories have a 'fun- 
damental' quasiparticle, which we call X, that acts as a vor- 
tex for the Z boson. X quasiparticles are non-Abelian anyons 
with quantum dimension \fm. We will call them metaplectic 
anyons, for reasons that we will explain. When two X parti- 
cles are fused, the result can either be the vacuum or one of 
a set of quasiparticles which we call Yj, with i = 1,2, ... ,r, 
and r = (m — l)/2. The Yi particles have quantum dimen- 
sion 2, but this does not mean that they are trivial; they are also 
non-Abelian anyons. Finally, there is a particle X' , which re- 
sults when X and Z are fused. Only a subset of the tensor 
product of the quasiparticles of the SO (171)2 TQFT and the 
quasiparticles /, a, ip of the Ising TQFT satisfy the constraint 
that they braid trivially with the electron = ip ■ Z, as we 
will describe in detail. We call the resulting topological phases 
metaplectic-Majorana TQFTs. 

A collection of quasiparticles of type X at fixed positions 
has an -dimensional degenerate state space in the SO(rn)2 
TQFT with tin ~ m N ^ 2 . Braiding these quasiparticles gen- 
erates unitary transformations in U (tin). These unitary trans- 
formations form a finite group, as in the case of Ising anyons, 
but unlike Fibonacci anyons. Therefore, it is not possible to 
make a universal quantum computer purely by braiding X par- 
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tides. We show that braiding can be efficiently simulated by a 
classical computer by showing that braiding operations satisfy 
a generalization of the Gottesman-Knill theore m 35 ' 36 . Indeed, 
the link invariants computed by these particles in a braiding 
process is known to be classically computable in polynomial 
time. However, the Yi particles - which one might naively ex- 
pect to be trivial since they have integer quantum dimensions 
- compute a link invariant (the Kauffman polynomial 37 at spe- 
cial points) that is #P-hard 3 -. Therefore, braiding Yi particles 
cannot be efficiently simulated classically. This does not mean 
that we can solve #P-hard problems since that would entail 
measuring the amplitude for a process with arbitrary accuracy. 
Indeed, as we show, the most straightforward approach to en- 
coding quantum information in Yi particles leads to a compu- 
tational model that can be efficiently simulated classically, and 
the image of the braid group of Yi particles is finite. Neverthe- 
less, the #P-hardness of braiding Yi particles hints that meta- 
plectic anyons and metaplectic-Majorana anyons may have 
computational power beyond a classical computer, in spite of 
the fact that they cannot serve as a universal quantum com- 
pute. In this respect, they may be similar to the linear optics 
model of Ref. [39l 

We will argue that our topological phase of metaplectic 
anyons is closely related to a set of recently proposed two- 
dimensional 40 and quasi-one-dimensional systems 4 ^ 4 -. In 
these systems, there are defects with interesting topological 
properties. In Ref. |4CJ they are dislocations in a fractional 
quantum Hall state in a Chern number 2 band. In Refs. I4ll - l43l 
the defects live at the edge of a fractional topological insulator 
or the edge between two v = 1/m quantum Hall states that 
are oppositely spin-polarized. There are two types of gapped 
edges, and a defect lives at the point-like boundary between 
the two types of gapped edges, generalizing the m = 1 case, 
in which they are Majorana zero modes. A form of braiding 
can be defined for the defects in these models. We show that 
this braiding operation is projectively equal to that of a ■ X 
quasiparticles in the metaplectic-Majorana TQFT However, 
there are important differences between metaplectic-Majorana 
anyons and the defects in these models, as we will discuss. 

We also note that related topological phases have been con- 
structed in Refs. I45l - l47i These topological phases have sim- 
ilar anyons with similar quantum dimensions and topological 
spins, but it is not clear what the precise relation is to our 
phases. 



n. SLAVE PARTICLE FORMULATIONS 

In this section, we give two slave particle descriptions of 
electronic systems in the topological phases that we discuss 
in the remainder of this paper. The first is a 'parton' model 4 ^ 
in which the electron operator is rewritten in terms of partons, 
each of which condenses in a simpler topological phase. The 
second is a non-Abelian analog of the flux attachment opera- 
tion that transforms electrons into 'composite bosons'^ 3 ^ or 
'composite fermions'— ~—. 

For later convenience, we fix notation for SO(m) represen- 
tations. We will often write m in the form m = 2r + 1. We 



use the standard notation that Ai, A2, . . . , A r are the funda- 
mental weights of SO(m). The representations with high- 
est weight Ai, A2, A3, . . . , A r _i, 2A r correspond to the rep- 
resentations of SO(m) on, respectively, vectors; two-index 
anti-symmetric tensors; three-index anti-symmetric tensors; 
. . . ; (r — l)-index anti-symmetric tensors; and r-index anti- 
symmetric tensors (with all indices running from 1 to m). The 
representation with highest weight 2Ai is the representation of 
SO(m) on two-index symmetric traceless tensors. The repre- 
sentation with highest weight A r is the spinor representation 
of SO(m). 

We first consider the following representation of the elec- 
tron annihilation operator: 

y c \x) = f(x)C aPX l(x)xl(x) (1) 

Here, /, Xa> an ^ Xa are fermions and a, /3 = 1, 2, . . . , 2 r . 
C a /3 is the intertwiner between two copies of the spinor rep- 
resentation of SO(m) and the trivial representation. This ex- 
pression for the electron is highly redundant, as is reflected in 
its U(l) x 0(m) gauge symmetry. The U(l) gauge transfor- 
mation is: 

f(x)^e™f(x), X 1 f(x)^e ie X 1 f(x) (2) 
while the 0(m) gauge transformation is: 

X 1 f{x)^O afj {x)xY{x) (3) 

We now suppose that the fermions / condense in a p + ip 
superconducting state while the fermions x^ 2 are m g a PP e d 
insulating states in which they fill a band with Chern num- 
ber equal to 1. Integrating out the fermions Xa' 2 > we generate 
a Chern-Simons term at level 2 for the SO(m) gauge field. 
(Note that we could, alternatively, consider a representation 
of the electron operator in which x a = Xa t> ut these fermions 
are in a gapped insulating state in which they fill a band with 
Chern number equal to 2.) Meanwhile, the excitations of a 
p + ip superconductor (coupled to a 2 + 1-D U(l) gauge field, 
which eliminates the Goldstone boson by the Anderson-Higgs 
mechanism) are those of the Ising TQFT. Naively, the excita- 
tions of this phase are simply those of SO(m)2 (which we 
will discuss in detail in the next section) tensored with those 
of the Ising TQFT. However, a vortex in the p + ip supercon- 
ductor of /-pairs will be accompanied with one half of a flux 
quantum in the Chern insulating states of x\' 2 ■ This flux will 
be produce a x\ 2 quasiparticle, carrying the spinor represen- 
tation of SO(m). Thus, a a quasiparticle in the Ising sector 
of the theory is accompanied by a quasiparticle in the spinor 
representation of SO(m). 

We now consider a (related and, possibly, dual) slave 
fermion description of an electron system in which we write 
the electron annihilation operator as: 

**(x)=f(x)z a (x) (4) 

where / is a neutral, spinless fermion and z a is a charge-e, 
spin- 1/2 boson, and a =t,|- 

We now rewrite the fields z a in terms of auxiliary fields 
in a non-Abelian analog of the flux attachment operation that 
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transforms electrons into 'composite bosons'^2i2i or 'compos- 
ite fermions'— ~—. This is simply a rewriting of the model, 
and the original model and the rewritten model would have 
the same solution if we could solve them exactly. However, 
this re-writing suggests an approximation that we might not 
otherwise consider. 

We replace the fields z a by auxiliary bosons Z a coupled to 
two SO(m)i Chern-Simons gauge fields, a 1 , a 2 . The fields 
Z a are m x m matrices that transform under SO(m) x SO(m) 
as Z^ — > C>2Z^O\ and Z^ — > O^Z^Oi, i.e. they transform in 
the fundamental representation of both SO(m)s. An SO{m)\ 
Chern-Simons gauge field would make Z a into a fermion. 
Therefore, two such gauge fields leave Z a a boson. In terms 
of these fields, the Lagrangian then takes the form 



C = Zi (id - a\ — ag) 



2mz 



(id - a] - a 2 ) Z t \ 2 



+ Z\ (id - a + a%) Z i + — !- I (id ~ a] + a 2 ) Zj 2 

+ p (id -ao)f+ ^-\(id + on) jf 
2m/ 

+ V(Z a ,f,p)+£ cs ( ai ) + £c S (a 2 ) (5) 

The relation between the original fields z a and the new fields 

Z a is: 



III. PARTICLE TYPES, TOPOLOGICAL SPINS, AND 
FUSION RULES 



We introduce the following notation for these quasipar- 
ticles. The particles carrying SO(m) representations A r 
and Ai + A,, will be called X and X' . The particles car- 
rying representations Ai, A2, . . . , A r _i, 2A r will be called 
Y\, Y2, . . . , Y r -i, Y r . Finally, the particle carrying 2Ai will 
be called Z. The particle carrying the trivial representation of 
SO(m) is equivalent to the vacuum from a topological point 
of view. We note that the special case m = 3 is equivalent 

to SU (2)4, and the X, Y\, X', Z particles correspond to spins 

1 1 3 o 

2 ' 2 ' 

The topological properties of the metaplectic TQFT are 
as follow a 49 ' 50 . The topological spins 6 a = e 2 ^ ha of these 
particles are given by hj = 0,hz = l,hx = %>hx' = 

r+4 



2m ' 



Their fusion rules are: 



X -X = I 



Y 



X ■ X' = Z + ]Tk, , 



X ■ z = 
X-Yi = 
Yi-Y = 



X', 

X + X' . 



Z-Y i =Y i 
Z-Z = l . 



j\ ~r~ ^nin(i+j,m— i— 3) 1 ^ 7^ 3 



K-K = I + Z + Y 



min(2i,m— 2i) 



(7) 



Zf(x) = Ve- l ^ a2 Z t (x)Pe i ^ ai 
Zi \ x ) = Ve lI ^ a2 Z l (x)Ve i ^ ai 



(6) 



We now assume that Z^ condenses, thereby breaking 
SO(m) x SO(m) to the diagonal SO(m). The Meissner ef- 
fect due to Z^ forces a* = a 2 , which we now write simply as 
a^. The two Chern-Simons terms then add, and a M has level 
2. 

We are now left with Zf, coupled to an SO{m)2 Chern- 
Simons gauge field. Decomposing Z^ into irreducible repre- 
sentations of SO(m), we have fields carrying the trivial rep- 
resentation, and the representations with highest weights A2 
and 2Ai. Since m(SO(m) x SO{m)/SO{m)) = Z 2 , there 
are also topological defects in the Z^ condensate. By forming 
combinations of the irreps in Z^ and the topological defects 
in Z±, we have particles carrying all of the allowed repre- 
sentations of SO(m)2, namely representations with highest 
weights 0, Ai, A2, . . . , X r , 2A r , Ai + A r , 2Ai. We will call the 
SO(ra)2 TQFT the metaplectic TQFT, for a reason to be ex- 
plained when we discuss quasiparticle braiding. 

The fermions / are assumed to condense in a p + ip paired 
state. Therefore, there are, in addition to the particles listed 
above, vortices a and fermions ip. This breaks the U(l) gauge 
symmetry / — > e l6 f,z — > e~' l8 z down to a Z2 symmetry. 
Consequently, a particles, which are vortices in the (//) con- 
densate are accompanied by Z2 flux which also inserts a topo- 
logical defect in the Z± condensate. As we will discuss in the 
next section, this means that only certain combinations of the 
particles in the Ising TQFT and the particles in the metaplectic 
TQFT are allowed. We dub this combination the metaplectic- 
Majorana TQFT. 



For the m = 3 case, there is a single Yi, which we will simply 
call Y = Y\, and the last of these fusion rules is modified to 
Y ■ Y = I + Z + Y . We obtain the dimensions of multi- 
particle Hilbert spaces from these fusion rules. If we denote 
the Hilbert space of n particles of type X with total charge Q 



X 



then 



<M?4 x) = hi™"' 1 ± 1) . A MuXn,x) 



,n-l 



<M^n£,x) = |(™ n ±l)- (8) 



Combining the Ising (see, e.g. Refs. l7lfITl) and metaplec- 
tic TQFTs, we naively have the particle types {/, a, ip} x 
{I, X, X', Yi, Z}. However, some of these are not local 
with respect to the electron operator = ip ■ Z. The 
topologically-distinct ones that are local with respect to the 
electron are: I,aX,ip,Yi, Z. These 4 + r particle types 
determine, for instance, the ground state degeneracy of the 
metaplectic-Majorana TQFT on the torus. However, it is 
worth noting that this is actually a Z 2 -graded TQFT, and 
one should also consider as distinct the particle types that 
differ from these 4 + r particle types by a single electron: 
ipZ, aX',Z, iPY,. 

Turning now to the particles allowed in the full metaplectic- 
Majorana TQFT, we have: 



(2m) 



n-l 



dim(H^5)-2"- 2 (m»±l) (9) 
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IV. F- AND i?-MATRICES 

We can determine the braiding properties of these particles 
using their F and R matrices. There are many non-trivial F- 
matrices for SO(m)2, which can be obtained by solving the 
pentagon identity. Some, which we will use below, are-: 



Consequently, we can represent the braid group in the fol- 
lowing way. We define the extra special group H(n,¥ m ) 
(sometimes called the Heisenberg group) generated by 
z, Mi, U2, • ■ ■ , u n , satisfying the relations 



F 



XYtYt 



X 



b X' 



pYiYiYi 




UiU i+ i 
UiUj 
UiZ 



1 , z m = 1 

ZU i+1 Ui 

j»t! \i — j| > 1 



UjU 

zu 



(15) 



(10) 



The F£ -matrix is an (r+1) x (r+1) matrix. For m = 3,5, 
it is given by, respectively^ 



This is a group of order m n+1 which is discussed further in 
Appendix lAl We introduce this group because, given a rep- 
resentation of H(n,¥ m ) by operators u, acting on a vector 
space, we can define a representation px of the braid group 
B n , as we we will see below and will discuss in further detail 
in Appendix lAl We construct a representation of H(n,¥ m ) 
of the requisite dimension as follows. Suppose, for the sake 
of concreteness, that n is even and that we are interested in 



ipXXX 

b x 



jpXXX 

b x 



1/V3 v/6 
V5 VTo 



'H I n . Then, we can define = span( fci,&2 



«n/2 



10 



hi € F m , and define the action of H(n, ¥ m ) on H^n by 
u 2l -i \ki,...,k n / 2 ) =u 2k ' \ki,...,k n / 2 ) 



)) with 



VlO -3(5 + V5) K 5 -^) P> U 2 i\h,...k n/2 )=\k 1 ,...,k i -l,k i+1 + l...,k n/2 ) 

) (16) 



10 2(5 -V5) -|(5 + V§) 



z \ki,. ..,ki,..., k n ) = 10 2 \ki, 



v n /2 



Similarly, the i?-matrices can be obtained by solving the 
hexagon identity. Some of the non-trivial ones, which we will 
use below, are: 



R 



xx 



:{_r-j)(r-j+l)-j „™(§+4 



We could have represented z by any 77i th -root of unity, but we 
have chosen uj~ 2 for later convenience. 

With this representation of H(n, ¥ m ) in hand, we define a 
representation px of the braid group B n according to: 



R YlYl = e 7ri ( m + 1 )/ m : ji 



R 



^7ri(m — l)/m &XZ 
5 ,Hvi 



Y-iYi _ TTi/m 



nX'Z 
K X 



X 1 



(12) 



1 



-(r 2 +r/2) ^ U fu{ 
j=0 



(17) 



With these F- and i?-matrices, we can compute how the 
states in the multi-quasiparticle Hilbert spaces of dimensions 
(O transform under braiding. 



V. iV-PARTICLE BRAID GROUP REPRESENTATIONS 

We now consider a situation in which we have n particles of 
type X in the SO(m)2 TQFT Braiding these particles leads 
to a representation px of the n-particle braid group, £>„. We 
now describe this representation and its image. Let px(o~i) be 
the representative of the braid group generator Ui (a counter- 
clockwise exchange of particles i and i + 1) acting on the 
n-particle Hilbert space. From the i?-matrices, we see that the 
eigenvalue equation for px{ai) is 



I (px(cTi) - fr-M'-i+V-i e™^ + ^> 
or, equivalently, 



(13) 



3=0 

where u = e 27ri / m . 



- r r2 oj j2 ) = 



(14) 



Direct computation shows that px (<7j) obeys the Yang-Baxter 
equation. Moreover, the states J^^Wi) are eigenvectors of 
the braid generator ( [17] ) with the same eigenvalues as Eq. ( TBI ) 
by virtue of the quadratic Gauss sum, uj^ ui^ = uj~ k ~ . 

The eigenvalues and dimensions determine the characters of 
the representation which, in turn, determine the representa- 
tion. Therefore, we conclude that ( fT71 ) is the representation 
( fT4b for n X-particles. This representation of the braid group 
is called the Gaussian representation^. 

We note in passing that there is another possible braid 
group representation on this Hilbert space, the Potts 
representation^, in which p(cri) = (t + 1)^ 2j=o u l ~ ^> 
and 2 + t + t^ 1 = m. The Potts and Gaussian representations 
coincide for m = 3, but differ for m > 5, where the Potts 
representation is not relevant to 50 (771)2 since the eigenval- 
ues of the braid group generators are different. Note that the 
m = 3 Potts representation is not related to the critical point 
of the ferromagnetic 3-state Potts model, which is the theory 
of Z 3 parafermions; it is, instead, related to the critical point 
of the anti-ferromagnetic 3-state Potts model^. 

The image of the braid group in the Gaussian representation 
can be understood as follows (see Appendix lAl and Refs. H^, 
I52I for further details). From Eqs. (TT3T > and ( TTTb . we see that 



[pxt^+i)] 1 ", px{o~i 



LU U i+1 Ui 



5 



[pxK-l)] 1 Ui Px(cTi-l) 



\i-j\>l (18) 



Therefore, braiding transforms any Uj into a product of ujs, 
up to factors of ui. If we mod out by the factors of u, then we 
have H{n — 1, m)/Z(H(n — 1, m)), the extra special group 
modulo its center. Braiding transformations are, therefore, au- 
tomorphisms of H(n—1, m)/Z(H(n—l, m)). Hence, the im- 
age of the braid group is a subgroup of the group of automor- 
phisms of H(n — 1, ¥ m )/Z(H(n — 1, m)). As we discuss in 
Appendix [A] this is equal to the metaplectic representation*^ 
of Sp(n— 1,F TO ) for n odd and Sp(n — 2,¥ m ) X H(n — 2, m) 
for n even. For this reason, we call X particles metaplectic 
anyons and we call SO(m)-2 the metaplectic TQFT. 

The group 5p(n — 2, F,„) x ii(n — 2, m) is a natural gen- 
eralization of the Clifford group. Recall that the Pauli group 
is composed of products of ± Pauli matrices for n/2 spins; in 
our notation, it is equal to H(n, 2). The group of automor- 
phisms of the Pauli group that are trivial on its center is the 
Clifford group, and it is equal to Sp(n, F 2 ) x V n /2- In other 
words, braiding metaplectic anyons generates a subgroup of 
the analogue of the Clifford group for qudits, with ¥2 — >• F m . 

Turning now to the full metaplectic-Majorana TQFT, we 
combine Eq. (fTTI i with the braid group representation for Ising 
anyons^ 

1 m— 1 

V * k=0 V „-_ n 



3=0 



(19) 

where vf = 1, ViVi+i = — Vi+\Vi, ViVj = VjVi for > 1. 



VI. QUANTUM INFORMATION PROCESSING WITH THE 
METAPLECTIC-MAJORANA TQFT 



the notation Yq = I) and, for r < j < m — 1, to the state with 
01 = V^m-j. a 2 = aX'. 

For such a qudit, there are two generators of the unitary 
transformations that can be performed by braiding. The first 
is a counter-clockwise exchange of the two crX-particles on 
the left. This implements the following gate which is diagonal 
in this basis: 



P(o-i) \j,n^) = e »i 



*W 3 |j,7ty) 



(20) 



The second is a counter-clockwise exchange of the middle 
two erX-particles. This can be obtained by using the F- 
matrix to transform into a basis in which these two particles 
have a fixed fusion channel, applying the R matrix, and trans- 
forming back, i.e. from F~ 1 RF. For the sake of concise- 
ness, let us consider the case m = 3. Then p{a^) \ j, n^) = 



M jk L 



k, n' b ) where 



|(1 + 2") 



M= I a^(i_ w ) i(2 + 




^(1-c) 0^ 



e s 
71 



1 -i 
-i 1 



(21) 

In a similar manner, we obtain the gate associated with a 
a counter-clockwise exchange of the last two aX -particles 
for m = 3, which takes the form p(a^) | j, n^) = 

Mj k e~^ 'e^ n ^ \k,n^) with 




I (1 + cj)/2 (l-w)/2 
(l-w)/2 (l + w)/2 ; 



(22) 



For multiple qudits, we can employ either a dense or sparse 
encoding. A dense encoding using 2k a X -particles can be 
represented by: 



aX aX aX aX 



aX aX 



We will consider three different encodings of quantum in- 
formation into the many-particle states of the Metaplectic- 
Majorana TQFT. For reasons that will become clear, we call 
them the 'qudit', 'qubit', and 'qutrit' encodings. 

Consider the state space of 4 aX-particles with total topo- 
logical charge Y\. It can be depicted graphically as follows. 



aX aX aX aX 



Y 



Oi 



a 2 



The first two particles fuse to 01, which can be 
I, Yi, . . . Y r , ijj, ipYi, ■ ■ ■ , i/jY r . In all of these cases, 02 = crX 
is possible. However, if a\ = Yi, . . . Y r , ipY\, . . . , ipY r , then 
02 = uX' is also possible. Therefore, there are 2(r + 1) + 
2r = 2m such states. We will take a basis \j,n^) with 
< j < m and = 0,1 for this 2m-state qudit. \j, 0) cor- 
responds, for < j < r, to the state with a\ = Yj, 02 = &X 
(with the notation Yq = I) and, for r < j < m — 1, to the state 
with ai = Y m _j, a,2 = crX' . Meanwhile, \ j, 1) corresponds, 
for < j < r, to the state with a\ = ipYj, 02 = oX (with 



Such an encoding uses 2k aX-particles for k— 1 qudits. How- 
ever, an exchange of neighboring particles will necessarily in- 
volve neighboring qubits. Consequently, simple single-qudit 
gates are complicated in terms of braids and errors in one qubit 
tend to infect others. We can, alternatively, use a sparse en- 
coding, such as: 



aX aX aX aX aX aX aX aX 



aX aX aX aX 



Yi 



In such an encoding, 4k aX -particles are used for k qubits. 
There are k sets of 4 aX -particles. Each set of 4 has total 
topological charge Y\ . These sets of 4 are paired so that each 
pair of sets (i.e. a group of eight aX -particles) has total topo- 
logical charge I. 

An alternative encoding scheme, which we call the qubit 
encoding, uses a aX particle and (n + 1) Y\ -particles (or any 
other Yi) to encode n qubits. It is depicted as follows: 
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Y x Yi Fi Fi Y x Y x 

aX 1 1 1 1 ■•• 1 1 aX 

a\ 0,2 03 0.4 a n 

where a; = aX or aX' . In order to express the gate that 
results when particles % and i + 1 are exchanged, it is useful to 
define = X, if m = 3 and Hi = Z^xX^Z^x if to > 5 
(note that Xi, Zi are the usual Pauli matrices here because 
we have qubits rather than qudits). We label the qubits by 
z = l, ...,n, and we define Zq = Z n+ x = +1. In addition, 
we define NOT~ + = / if Z l - 1 Z i+1 = 1 and NOT~ + = Xi 
if = —1. Then, a counter-clockwise exchange of 

particles i and i + 1 results in a gate that can be written in the 
following form: 

py I (<7 i ) = e£ Hi NOTy+ (23) 

Finally, we introduce one more encoding: the qutrit rep- 
resentation. (Qutrits are obtained for any to. Note that the 
qudit representation introduced earlier is never a qutrit repre- 
sentation since the dimension 2m is always even). A qutrit is 
encoded in four Y\ particles with total charge /: 

Yx Y x Y x Y x 

I 1 1 1 1 / 

a 

From the fusion rules ©, we see that the charge a = I, Y 2 , Z 
(except in the case to = 3, where a = I,Yx, Z). Braiding the 
first two particles enacts the transformation: 

(l 0\ 

p Yl (ax) = -e™' m ZJ (24) 

v° -v 

while braiding the second two enacts: 

/ 2uJ 2^2 -25J\ 
PyM = [2V2 2\/2 (25) 
\-2uJ 2V2 2uJ J 



VII. CLASSICAL SIMULATION OF BRAIDING IN THE 
METAPLECTIC-MAJORANA TQFT 

Regardless of the encoding, universal quantum compu- 
tation is not possible purely through braiding because the 
braid group representation dTTb for n A -particles is contained 
within Sp(n — 1, F m ) for 11 odd and Sp(n — 2, F m ) k H (n — 
2, ¥ m ) for n even. As we discuss in greater detail in Appendix 

n 

\S P {2n, ¥ m )\ = to" 2 H (to 24 - 1) (26) 

i=X 

while \H(n, m)\ = m n+1 , so the braid group has a finite im- 
age under the Gaussian representation. Therefore, it is not 
possible to approximate an arbitrary unitary transformation to 



any desired accuracy. In fact, braiding aX particles can be 
efficiently simulated by a classical computer. 

Since it is known that braiding in the Ising TQFT can be 
efficiently simulated classically, we focus on the braiding of 
metaplectic anyons. Recall that braiding metaplectic anyons 
transforms products of w^s into products of UiS, as we noted 
in Section QVj As a result, the evolution of eigenstates of such 
products can be efficiently simulated classically by follow- 
ing the evolution of these operators. In order to see this in 
greater detail, it is convenient to embed H (n — 2, F m ) inside 
H{2n, F m ) as follows. Let Xx, ■ ■ ■ , X n , Zx, ■ ■ ■ , Z n , to be a 
set of generators of H(2n,¥ m ), as described in Appendix lAl 
(see, especially, Eq. IA3b . Then Ui — XiXi + iZiZ\ +1 faith- 
fully represents the extra special group (fT~5b . Consequently, 
Px(cj) = -^7= Sj'l~o ^ U$ represents the braid group. We 
can prepare states that are eigenstates of Ui by creating pairs 
out of the vacuum. Such states are stabilized by products 
of Xi and Zj operators since Ui can be expressed as such 
a product. To see how any state stabilized by products of 
Xi and Zj operators transforms under braiding, we can fol- 
low the evolution of the operators Xi, Zj. It is sufficient to 
consider the case of two qudits. We would like to see how 
Xx, X2, Zx, Z2 (and, therefore, the group that they generate) 
evolve under the action of R. First, note that we can replace 

the set Xx,X 2 ,Zx,Z 2 by the set Zx, XxX 2 , Z X Z\, X x Zx, 
which generates the same group. The latter three commute 
with U and, therefore, with R. Therefore, we need only 
study how Zx evolves. Using W Zx = u~^ZiU 3 ', we see 
by direct computation that px{ux)Zxp x {<7x) = u~ k ZxU k , 
where k = (to + 1) /2. Therefore, the evolution of Z\ can be 
efficiently simulated classically and, as a consequence, so can 
the evolution of any state stabilized by products of Xi and Zj 
operators. Thus, we conclude that we can efficiently simulate 
classically any operation that consists of creating pairs of X 
particles out of the vacuum, braiding them, and then measur- 
ing them a basis of products of X{ and Zj operators (e.g. the 
U 2 i-x basis). 

Of course, as noted above, H(2n,¥ m ) is much too 
large. It associates an TO-state qudit to each AT -particle 
while, in the dense encoding, there should be a qudit as- 
sociated to each pair of X-particles. Therefore, braiding 
should commute trivially with roughly half of the gener- 
ators of H(2n,¥ m ). This is, indeed, the case, as may 
be seen by considering the following set of generators 
of H(2n, F m ): U u ... U n -x,Ux, . . . U n -x,XxZx,X n Zl,w 
where {/, = XiXi + xZ^ Z- l+ x - The generators C/^s, XxZx, and 
X n Z\ all commute with the L/jS and, therefore, with braiding. 

Braiding is not universal in the qubit representations, either. 
We now show that the group generated by the py ± {<Ji) oper- 
ators acting on the qubit representation is finite, and we give 
an efficient classical way to store an arbitrary element of this 
group and to efficiently compute products of elements of this 
group with braid generators (the method we describe will only 
store an element up to an overall phase). For all to (including 



7 



both m = 3 and m > 3), a direct computation gives 

(NOT7'+) HiNOTr^ = (27) 
(nOT7/' + ) iT; +1 NOT7/' + = ffcfli+i, 
(N0TT'+) t i? i _ 1 N0T7' + = ff^JJi, 

so conjugating a product of the Hi by a unitary NOT" ,+ gives 
some, possibly different, product of the Hi. The group gen- 
erated by the operators e™ Hi is an Abelian group, which we 
call G. Since e 2lrHl = 1, we can write an arbitrary element 
of the group as e 1 ^** ki ik H \ where the ki are integers ranging 
from 0, 2m — 1, so the group is a subgroup of IQm- How- 
ever, since e lrHi = — 1, there are only 2 • m™ distinct group 
elements which can be written as (±1) • e l ^ i ki ™ H \ where 
the ki are integers ranging from 0, m — 1. This group is in 
fact Z™ 4 x Z 2 , and the generators of the group can be taken 
to be —e % ™ Hi and —1. The group generated by the operators 
NOT~' + is a subgroup of the Clifford group; call this group 
H. Then, because conjugation by NOT~' + defines an auto- 
morphism of G, the group generated by e l ™ Hi NOT~ + is the 
semi-direct product G x H. This gives us an efficient way 
to store elements of the group by storing a list of integers fc; 
and also storing an element of the Clifford group. We spec- 
ify an element U of the Clifford group by specifying UXiW 
and UZiW for all i. These products UX^ and UZ { W are 
products of Pauli matrices and so can be stored efficiently (we 
are essentially using the Gottesman-Knill theorem here). Stor- 
ing these products fully specifies UOU^ for any operator O 
and so specifies U up to a phase. To take a product of two ele- 
ments of the group, say the first being represented by a product 
AU and the second by a product A'U' where A, A' are in the 
Abelian group and U, U' are in the Clifford group, we write 
AUA'U' = A(UA'W)UU'. We then compute UA'W using 
our known values of UXiW and UZiW and the result will 
be some other element of the Abelian group, all it A", Then 
the desired product is AA"UU', and the product of the first 
two is in the Abelian group and the product of the second two 
is in the Clifford group. 

It should not be surprising that the group image is finite. 
The Y\ particles can be obtained by fusing a pair of X par- 
ticles. Thus, the fusion tree in Section [VT] that defined the 
qubit representation can be written as a tree with 2(n + 1) + 2 
AT -particles, with 2(n + 1) of the .A -particles fusing in pairs 
to make (n + 1) Y\ -particles. Braiding two Y\ -particles can 
be done by braiding two pairs of A -particles. Since the im- 
age for braiding A -particles is finite, it is no surprise that the 
image for braiding Y\ -particles is also finite. However, it is 
still important to check, as we have done, that we can effi- 
ciently store elements of this group; after all, the tree that we 
have written here with 2(n + 1) + 2 A -particles is related by 
some sequence of F moves to the previous tree in terms of 
A -particles and it is not immediately obvious that all these F 
moves can be computed efficiently. 



VIII. COMPUTATIONAL COMPLEXITY OF LINK 
INVARIANTS 

In the previous section, we have seen that braiding is 
not universal for quantum computation in any representation. 
Moreover, braiding in the qudit and qubit representations can 
be efficiently simulated classically. However, this theory dis- 
plays a surprise when we turn to the computation of link in- 
variants. Thus far, the most-studied examples of TQFTs for 
which braiding is universal for quantum computing have been 
precisely those for which an evaluation of the link invariants 
is $P-hard. However, there seems to be no deep reason 
why this should be true generally, and indeed the present the- 
ory is not universal for quantum computing (through braiding 
alone), but it does have a link invariant that is ^P-hard to 
compute. Said differently, there are experiments whose re- 
sults are #P-hard to predict, i.e. cannot be predicted with a 
classical computer (unless the hierarchy of complexity classes 
collapses), even though braiding alone is not sufficient for uni- 
versal quantum computation. 

We give a more precise definition of this link invariant 
elsewhere^. Here we will give its physical motivation. We 
imagine creating a collection of pairs of Y\ particles out of 
the vacuum. We braid them with each other and then fuse 
them again in pairs. There will be some amplitude E(L) for 
all of these fusion processes to give the vacuum, i.e. to be 
annihilation processes. (When two Y\ particles are fused, the 
result could be the vacuum /, but it could, instead, be Z or 
Y2, except in the m = 3 case, in which there is no Y% particle 
and it could, instead, be Y\ .) Here, L is the link formed by 
the spacetime trajectories of the Y\ particles. The amplitude 
E(L) is our 'link invariant'. We use quotation marks because 
this amplitude is not necessarily a topological invariant unless 
further conditions are satisfied. However, if the interaction be- 
tween the Yi particles decays exponentially (or faster), then, 
in the limit that the particles stay far apart while braiding, this 
amplitude will depend only on the topological class of the Y\ 
trajectories. When the particles are being pair-created and an- 
nihilated, the amplitude will acquire a non-topological, non- 
universal phase. However, this can be made to cancel between 
creation and annihilation. Alternatively, if two different braid- 
ing processes are interfered, then this non-topological phase 
will cancel. See, for instance, Refs. 1X111561 for a discussion of 
interference measurements for link invariants. 

The starting point for the #P-hardness of E(L) is a result 
of Lickorish and Millett^i. They show that the link invariant 
E(L) can be written as 

E(L)= 5>- 4 < s < L - s \ (28) 

SCL 

where 

a = — i exp(— in/m). (29) 

Here, the sum is over links S which are a sublink of link L. 
A link may be made of more than one disconnected compo- 
nent, where each component of the link is some knot; we use 
c(L) to write the number of components of L. A sublink S 
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contains some subset of the components, so there is a total of 
2 C ( L ) terms in the sum, with each factor of 2 coming from the 
choice of whether a given component is in a sublink or not. 
We can specify a sublink S by a vector s with entries for 
i = 1, c(L), such that s; = +1 if the i-th component is in 
S and Si = — 1 otherwise. The invariant (S, L — S) is defined 
to be the sum of over pairs i E S and j e L — S, where 
is the linking number between the i-th sublink and the 
j-th sublink^ 

Eq. d28| i looks very much like the partition function of an 
Ising model at an imaginary temperature. The sum over sub- 
links corresponds to a sum over the "Ising spin" degrees of 
freedom s,;, while the term a~ 4 ^ S ' L ~' s ^ looks like a complex 
Boltzmann weight. To see this, write 

-4{S,L-S) = -£(i + -,,.)<*, j) 

= -2^(l-s iSj )(U) (30) 



Consequently Eq. d28l ) is equal to 

E{L) = a - 2 £ !<3 <^> a 2 ^<i s ^ {i ' j) . (31) 

se{-i,i} c < t ) 

So, up to the prefactor in front, the resulting link invariant is 
the partition function of an Ising spin system with Boltzmann 
weights 



(32) 



where j3 = —2iri/m + iri. 

Note that the temperature is purely imaginary. The quantity 
(i, j) plays the role of a matrix of coupling constants; note that 
these linking numbers can be taken to have any integer 
values. In particular, the Ising model need not be planar, and 
any choice of can be realized by some link L in which 
the number of crossings is at most polynomial in I {h 

We will now show that there is a class of links for which we 
can relate this Ising model with complex Boltzmann weights 
to more familiar models with real or even real and positive 
Boltzmann weights. We then argue that computing the result- 
ing partition function is #P-hard. (While we cannot relate 
E(L) for an arbitrary link to an Ising model with real Boltz- 
mann weights, it is sufficient to do so for the class of links dis- 
cussed below. We can then conclude that if we can compute 
E(L) for an arbitrary link, then we can solve any problem in 

To obtain an Ising model with real or even real and pos- 
itive Boltzmann weights, we use the following trick. We 
consider links L constructed as follows. We begin with a 
link L' with c(L') = N unlinked components, i.e. for any 
i,j G {1, 2, . . . , N}, = 0. We then add components 

N + 1, N + 2, . . . , c(L) to form the link L. They are cho- 
sen so that if i,j G {1,2, ... ,N}, then (i, k) = (j,k) (if 
k G {1,2,..., N}, then both sides of the equality are zero, 
but if k G {N+ 1, N+ 2, . . . , c(L)}, then they might be non- 
zero). We now evaluate the link invariant E(L) in two steps. 



First, we sum over the choices of for k = N +1, . . . , c(L) 
to define an "effective Boltzmann weight" for the first N Ising 
spin variables. Summing over component k generates an ef- 
fective interaction between i and j if (i, k) = (j, k) ^ 0. The 
effective Boltzmann weight will be real and E(L) is equal to 
the sum over the 2 N choices of the first N spin variables using 
the effective Boltzmann weigh 

Consider a pair i, j with 1 < i < j < N. We now add 
a component k G {N + 1, N + 2, . . .}, such that (i, k) = 
(i> k) = d for some d and such that (k, I) = for I different 
from i or j. Then, summing over Sfe = ±1 will produce an 
effective interaction between s, and Sj . Summing over Sk = 
±1 gives a weight 



^ a 2s,s k (i,k)+2s j s k (3,k) = 



2d(si-\-Sj )sk 



s k e{-i,i} 



where 



s fe e{-i,i} 

{Vy) s%S3 V~z, (33) 



and 



a -id + a Ad 



z = 2{a- id + a id ) 



(34) 



(35) 



and any ambiguity in the sign of the square-root is resolved 
by choosing -Jy^fz = aT Ad + a +Ad . 

Ignoring the overall factor y/z, the effective weight is 
(y/y) SiSj ■ By adding additional components k, k' of the link 
and summing over Sk>,Sk", and so on, we can replace this 
weight with any power, so that the effective weight for the 
first TV variables can be chosen to be (again up to an overall 
factor) 



l<i<j*<JV 



(36) 



for any matrix with non-negative integer entries (in fact, it 
is also possible to obtain negative entries by a slightly different 
trick but we will not need that here). The size of the link 
needed to produce this effective weight is at most polynomial 

The quantities y are real. However, depending upon m 
and d, they may be positive or negative. In fact, for any odd 
m > 1, we can choose — 1 < y < by an appropriate choice 
of d, and for odd m > 3 we can instead choose < y < 1 
by an appropriate choice of d. One way to obtain positive 
weights for m = 3 is to pick the entries of Jy to be even 
integers. In this way, we succeed in constructing a link in- 
variant that equals, up to multiplication by a trivial overall 
constant, the partition function of an Ising model at real, posi- 
tive temperature with antiferromagnetic couplings. By taking 
these couplings large, we can ensure that ground states pro- 
vide the dominant contribution to the partition function. That 
is, that the partition function is equal to No exp(—{3Eo) plus a 
small correction (small compared to exp(— /3Eq)), where (3 is 
now real and positive and where Eq is the ground state energy 
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and Nq is the number of ground states. Making the correc- 
tion small compared to exp(— j3Eo) requires only polyomi- 
ally large coupling constants (we are choosing the coupling 
constants large enough that energy outweighs entropy and so 
the sum of the weights of all the higher energy states is small 
compared to the weight of a single ground state). Then, an 
evaluation of the partition function lets one determine both 
the ground state energy and also the number of ground states. 
Counting the number of ground states is equivalent to finding 
the number of maximum cuts in a graph which is a #P-hard 
problems. Indeed the definition of #P is that it is the prob- 
lem of counting the number of solutions to a decision problem 
iniVP. 

This approach shows that evaluation of the link invariant to 
exponential accuracy is $P-hard. In fact, it is possible also to 
consider the case with negative and real Boltzmann weights 
(the case y < but has odd entries). Then, even the 
evaluation of the sign of the partition function is #P-hard, 
as follows from a result of Goldberg and Jerrum^. The sign 
of the partition function is equal to the phase of the link invari- 
ant multiplied by some overall phase which can be computed 
trivially. 

Similar behavior is seen in the theory ofS, which also has 
a finite braid group image but #P-complete link invariants. 
It would be interesting to see if our theory follows the pattern 
of their theory, where different approximations to the link in- 
variant are in P, or are SBP-hard, or are #P-hard, depending 
upon the accuracy of the approximation. It would be interest- 
ing to see if their theory, like ours, is classically simulable for 
certain measurements 



IX. RELATION TO FRACTIONAL QUANTUM HALL 
DEVICES 

In Refs. I4II - I43I (see also Ref. |44b , a model was presented 
in which the boundary of a fractional topological insulator 
was divided into 2N intervals, with the i th interval lying be- 
tween points Xi-i and Xi and xq = x^n- The even intervals 
(x2j-i, X2j) are brought into contact with s-wave supercon- 
ductors, while the odd intervals (x2j, %2j+i) are brought into 
contact with ferromagnets. The points x% are viewed as parti- 
cles. They 'fuse' to the 2m possible allowed total spins (mod- 
ulo 1) on the even intervals or 2m possible allowed charges 
(module 2e) on the odd intervals. They can be 'braided'^^ 2 - 
by a measurement-only proces o 60 ' 61 . The resulting unitary 
transformation for braiding two neighboring defects at Xk, 

U k , k+1 = e l ^' 2m (37) 

where q = 0, 1, ... , 2m — 1 are the possible charges/spins on 
the interval between the two defects. If we write q = mqj + 
2jM where qi = 0, 1 and jnj = 0, 1, . . . , m — 1, then 4 ^ 

U k , k+l = e" m ^ w& (38) 

where w = e 27 ™/"\ The first factor is the braiding transforma- 
tion for Ising anyons if m = 1 (mod 4) and for the opposite- 



chirality version of Ising anyons if m = 3 (mod 4). The sec- 
ond factor can be rewritten using the Gauss quadratic sum as: 

<J™ = -= V u?\{ (39) 
v j=o 

which is the same, up to a phase, as Eq. ( fTTb (see also Eq. 25 
of Ref. SI). 

Thus, these physical models give a very natural interpreta- 
tion to the elements of the extra special group H(2n — 1, F m ): 
these are the operators that rotate the phase of the super- 
conducting order parameter or the ferromagnetic spin by 
4/T. Their eigenvalues are just the allowed charges/spins on 
gapped intervals modulo charge 2e or spin-1. 

However, it is also important to note the differences be- 
tween the metaplectic-Majorana TQFT and the models of 
Refs. I4II - I44I The latter models are gapless since they 
have the Goldstone boson associated with superconductiv- 
ity (which is not given a gap by the coupling to a 3D elec- 
tromagnetic field). Therefore, these models are, at best, in 
quasi-topological phases^ and are related to the metaplectic 
TQFT in the same way that chiral p-wave superconductors 
are related to Ising anyons: they have some but not all of the 
properties of a true topological phase. Furthermore, we note 
that the models of Refs. I41M44] do not appear to have a Z 
particle. They have 2n-particle Hilbert spaces of dimension 
(2m)™ _1 . This is the same as the direct sum %2n ® W-2n< 
which suggests that these models do not distinguish between 
the Z particle and the vacuum. Moreover, the Yi particles 
are non-Abelian in the metaplectic-Majorana TQFT, but the 
charges/spins are Abelian anyons in the models of Refs. 
H . In the metaplectic-Majorana TQFT, when a Yj particle 
is taken around a Y k particle, a phase e ±m J k / m results, de- 
pending on whether they fuse to Y\j_ k \ or Y min ( i+ j m _ i _jy 
Each of these possibilities occurs twice (for each pair) if we 
allow the total charge to be I or Z. In the models of Refs. 
I411 - I441 however, the phase e m i k / m results when a charge j is 
taken around charge k or m — j is taken around m — k while 
e -mjk/m msu \i s wn en a charge j is taken around charge m—k 
or m — j is taken around k. In our model, we can only deter- 
mine the phase resulting from a braid by performing a mea- 
surement of the total topo log ical charge of the two particles. 
In the models of Refs. I4ll444l however, we can determine the 
phase resulting from a braid by simultaneously measuring the 
charges of the two intervals. 

A possible path to understanding the relation between the 
metaplectic-Majorana TQFT and the models of Refs. I4ll - l44l 
is through Slingerland and Bais' 6 - analysis of SU (2)4, which 
is equivalent to 5*0(3)2. They show that the condensation of 
the spin-2 particle (the Z particle), causes the confinement of 
the spin-1/2 and 3/2 particles (the X and X' particles). The 
Yi particle splits into 2 particles which, together with /, form 
a Z3 multiplet. A version of this scenario should occur for 
general SO{m)2, and may be related to the models of Refs. 
I4ll - l44t the charges/spins on intervals are the Abelian quasi- 
particles of the theory, which are the only 'true' quasiparticles 
in the theory since they are not confined, while X particles are 
confined but, if the energy required to pull them apart is sup- 
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plied, then a projective remnant of their non-Abelian braiding 
properties survives. The dislocations of Ref. [4oK63l may have 
a similar relation to the X particles of the metaplectic TQFT. 



X. DISCUSSION 

It was recently realized that the transformations associ- 
ated with Ising anyons could also be realized in three spa- 
tial dimension o 64 ! 65 . Although there is no braiding in three 
dimensions, extended objects, which could be viewed as par- 
ticles connected to ribbons, would have the topology of their 
configuration space governed by an enhancement of the per- 
mutation group, E(Z'i, n ~ 1 xi £271) (here, the E{...) denotes the 
restriction to elements whose combined parity is even). The 
Z 2 factors keep track of the twisting of the ribbons, modulo a 
Att twist, which can be undone. Solitons supporting Majorana 
zero modes realize a projective representation of this group, 
which has image H(n — 2, F2) XI S 2 n- Thus, the non-Abelian 
statistics of Ising anyons can be understood as simply permu- 
tations together with 2ir ribbon twists of pairs of particles. 
Two such twists anti-commute if they share a particle (but not 
both). The non-Abelian statistics of X particles in SO(rn)2 is 
a generalization of this to fractional twists: H(n—2, ¥2) is re- 
placed by H(n — 2, F m ) so that the (purely fictitious) ribbons 
connecting particles can be twisted up to m — 1 times. 

Although the resulting unitary transformations are richer 
than those of Ising anyons, this TQFT is still incapable of 
performing universal quantum computation through braiding 
alone. The braid group has an image which is finite. How- 
ever, a certain link invariant associated with the amplitude for 
creating pairs of Yj. particles, braiding them, and annihilating 
them in pairs is ^P-hard to compute. This suggests that there 
may be greater computational power lurking just beneath the 
surface of this theory and, perhaps, that it becomes apparent 
when braiding is supplemented by measurement at interme- 
diate steps of a computation. Specific protocols by which 
universal quantum computation could be achieved with meta- 
plectic anyons (with or without Majorana zero modes) are an 
interesting open problem. 
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Appendix A: Sp(2n,¥ m ) and the braid group of metaplectic 
anyons 

In this appendix, we will discuss in greater detail the im- 
age of the representation of the braid group associated with 



X -particles. We begin with a 2n-dimensional vector space 
V^n over F TO (with m assumed to be prime) equipped with 
a non-degenerate symplectic form [, ]. We can take as a ba- 
sis of this vector space Vj = (0, . . . , 0, 1, . . . , 0), which 
has zero for every entry except for the i th , which is 1. We 
will take the symplectic form to be [vi, Vj] = ±5i±ij. The 
group of linear transformations that preserve the symplectic 
form [, ] is the symplectic group Sp(2n, F m ). This is a finite 
group whose order can be determined as follows. We want 
all ways of choosing vi, . . . , i>2 n so that [v%> Vj] = ±<5i±i,j- 
There are m 2n vectors in V2 n since it is composed of all linear 
combinations of v\, . . . , i>2 n with coefficients in F TO . There- 
fore, there are to 2 ™ — 1 ways to choose v\ 7^ 0. There is 
a (2n — 1) -dimensional space of vectors v with [vi, v] = 0. 
Therefore, there are to 2 ™ — to 2 ™ -1 choices of vector V2 with 
[f 1 1 v 2\ 7^ 0. Since the possible non-zero values of [v\ , V2] are 
1, 2, . . . , m - 1, there are (to 2 ™ - to 2 "~ 1 )/(to - 1) = to 2 ™" 1 
choices of V2 with [vi, V2] = 1. Continuing in this way, we 
find that there are 

n n 

J (to 2 ' - 1) to 21 - 1 = m™ 2 [J (to 21 - l) (Al) 

i=l i=l 

elements of the group Sp(2n, F m ). 

Now consider Van as an additive group. Consider a cen- 
tral extension G: 1 — s- F m — > G — >• V2 n —> 1. Since 
V2n is Abelian, the commutator map G x G —¥ G given by 
92) —> <7i32<?i _1 <?2 _1 takes values in the center F m and is 
unaffected by multiplication by the center, so it defines a map 
V2n x V2n — > F m . In the case of the specific central extension 
that is usually called the 'extra special group' or 'Heisenberg 
group', which we denote by H(2n, F m ), this map is just the 
symplectic form [, ]. The elements of H(2n, F m ) can be writ- 
ten in the form (v,k), where v € V2n an d k E F m . The 
multiplication rule is (v\ , hi) ■ (ui , k\ ) = + V2, k\ + + 
[vi, ^2])- For the basis taken above with [v^, Vj] = ±§i±i t j, if 
write Ui = (vi, 0) and z = (0, 1), then we have the defining 
relations introduced in Sec. [V] 



UiUi+\ = zu i+ iiit 

UiUj = UjUi , \i — j\ > 1 
UiZ = zui . (A2) 

If we, instead, take a basis fi of Vm with [/2i-i,/2j] = 
Sij and [/ 2 i-i,/2j-i] = [f2i,hj\ = 0' then we have a 
different generating set for H(2n,¥ m ): Xi = (/ 2 ,_i,0), 
Zi = (f 2 i, 0), z = (0, 1) satisfying: 

XiXj = XjXi , ZiZj = ZjZi 
XiZj = z^ZjXi 
XiZ = zXi , Z t z = zZi (A3) 

These two presentations of H(2n, F m ) are related by u 2 i-i = 
Xi, u 2 i = ZiZ\ +1 for i 7^ n and u 2n = Z n . 

The symplectic group Sp(2n,¥ m ) of T^n acts on 
H{2n, F m ) in the natural way. These are automorphisms that 
act trivially on the center Z(H(2n,¥ m )) of H(2n,¥ m ). In 
addition, the inner automorphisms - conjugation by elements 
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of H(2n,¥ m ) - are also trivial on Z(H(2n,¥ m )). In fact, 
the group of automorphisms of H(2n, F m ) that are trivial on 
Z{H{2n,¥ m )) is given by Sp(2n,¥ m ) k V 2n . (V 2n , rather 
than H(2n, F m ), is the second factor in this semi-direct prod- 
uct because Z(H(2n,¥ m )) acts trivially on H(2n,¥ m ) by 
conjugation, so only H(2n,¥ m )/Z(H(2n,¥ m )) = V 2n ap- 
pears). The group Sp(2n, F m ) x H(2n, ¥ m ) is, therefore, an 
extension of the group of automorphisms of H(2n, F m ) that 
are trivial on Z(H(2n, F m )); the group has been extended by 
Z(H{2n,¥ m )). 

This is a useful extension to consider because, given an ir- 
reducible representation M of H(2n,¥ m ), there is a unique 
induced representation X of Sp(2n, ¥ m ) x H(2n 1 F m ) whose 
restriction to H(2n, F m ) is M, as shown in Ref . p4l and as we 
discuss in the next paragraph. Moreover, given a representa- 
tion AO) = u k of Z(if(2n,F m )), there is a unique induced 
representation A/ of H(2n, F m ) whose restriction to its center 
is A(/c). Here, ui is an m th root of unity. Let M v = M(v, 0) 
for (v,0) G H(2n,¥ m ). Then, the induced representation of 
H(2n,¥ m ) must satisfy M U M V = X([u,v])M u+v . Conse- 
quently, M™ = 1, M Vz M Vz+1 = u~ 2 M Vi+1 M Vi , M Vi M Vj = 
M Vj M Vt for \i- j\ > 1. 

This representation of H(2n, F m ) induces a representation 
of Sp(2n, F m ) k H(2n, F m ) as follows. Consider the action 
of # G Sp(2n, ¥ m ) on h £ H(2n, ¥ m ) by conjugation inside 
5p(2n,F m ) k ff(2n,F m ): ft, -> p/ig" 1 . Since #(2n,F TO ) 
is a normal subgroup of S'p(2n,F rn ) k iJ(2n,F m ), ghg^ 1 6 
H(2n,¥ m ). Therefore, for each g € Sp(2n,¥ m ) there is a 



representation of H(2n,¥ m ) given by ft, — > M, 



3^9 



But 



since there is a unique representation, there must be a unitary 
transformation X(g) such that M ghg -i = X(g)MhX(g)~ 1 . 
This defines up to a scalar. In fact, X(g) is not quite 

a linear representation of Sp(2n, ¥ m ). It is a projective rep- 
resentation or, equivalently, it is a linear representation of the 
double-cover of Sp(2n,¥ m ), namely the metaplectic group. 
This representation can be given explicitly in terms of the M v 
according to the relation: 



X{g)= Y, a v (g)M v 



(A4) 



where V\(g) = im(l — g). It may further be shown that 
0-vi.g) = X([u,g(u)]) a (g) where v = u- g(u) G Vi(s). 

We now consider the following map^i from B 2n+ i — > 
Sp(2n, ¥ m ). To the generator <7j of B2 n +i> we associate the 
Sp(2n, ¥ m ) transformation cr, that acts on V 2n according to 



<Ti(Vi) 



Vi 

Vi±i T Vi 

Vj, 



j\ >1 



(A5) 



It may be directly checked that this transformation preserves 
the symplectic form [, ] and that &i satisfy the defining rela- 
tions of the braid group. Then, from Eq. IA4I there is a braid 
group representation 



X(cr t ) = Y av(<n)M v 



(A6) 



From Eq. (IA5b . we see that Vi(cj) = {kvi\k e F m }. 
Hence, for g = cr^, we have v = kvi e Vx(<7,) and 
u = such that v = u — g(u). Consequently, a v (&i) = 

X([u,ai(u)})a (ai) = X([kv i+ i,-kvi])a ((ri) = w fc "a (cr i ) 
Therefore, 



M, 



kv. 



m—1 



k=0 



k=0 



Ml 



(A7) 



where M is a normalization constant. We see that this is the 
same as the braid group representation in Eq. ( ITTb which de- 
termines the braiding of 2n + 1 X-particles. Therefore, the 
image of the braid group representation of 2n + 1 AT -particles 
is equal to the metaplectic representation of Sp(2n, ¥ m ). 

The technical reason why the case of 2n + 1 particles is 
simple is that the braid group -Ban+i nas an even number 
of generators <ti , . . . , a 2n (since m exchanges particles % and 
i + 1). For an even number of generators, there is a natural 
mapping to Sp(2n, ¥ m ) since the latter is defined on a sym- 
plectic vector space, which must be even-dimensional. For an 
even number 2n of particles, the braid group B 2n has an odd 
number of generators. In order to construct the correspond- 
ing symplectic group, we begin with the symplectic vector 
space V 2n over F m and pick a vector e\ G V 2n . Then we 
consider the group G of linear transformations that preserve 
the symplectic structure [, ] on V 2n and leave e\ invariant. 
The vector space orthogonal to e\ is (2n — 1) -dimensional, 
so G is the odd-dimensional analogue of a symplectic group 
and is sometimes called an odd symplectic groups. Clearly, 
Sp(2n — 2, F m ) C G. The rest of G is given by transforma- 
tions of the following form. Let e 2n be the vector that satis- 
fies [ei,e 2n ] = L Then, for any v G span(e 2 , . . . , e 2n -i) 
and k G F m , the symplectic form [, ] and e\ are left in- 
variant by the transformations e 2n — > e 2n + v + ke± and 
ei — > Ci+ [v, ej]ei for i = 2, 3, . . . , 2n — 1. These transforma- 
tions, parametrized by (v, k) form the group H(2n — 2, m), 
as discussed above. They can be written explicitly in matrix 
form as 



/ e 2n \ 



G 2 





-f>n-l c 





6„_i 





In-l 0>n-l 


Vo o 


1 



f e 2n ^ 



\ ei I 



(A8) 



where a n -\, 6 n -i are ( n ~ 1) -component column vectors over 
F m , c G F m , I n -i is the (to— 1) x (n— 1) identity matrix, and 
the basis e 2 , e$, . . . e 2n -\ is chosen so that [e,-, e 2n +\-i\ = 1 
for i < ri and [ei,ej] — for j ^ 2n + 1 — i. This is 
precisely the group H(2n—2, to) in its representation as upper 
triangular matrices. Then, following the steps given above 
for an odd number of particles, we obtain a mapping B 2n — > 
Sp(2n - 2, F m ) x H(2n - 2, to). 
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There are two versions of these theories which differ by the 
Froebenius-Schur indicator, which accounts for the minus sign in 
these f-matrices as well as a few other differences. 
The linking number of two closed oriented curves may be com- 
puted by drawing a projection of the link on a plane, with care 
taken to denote over-crossings and under-crossings. Those cross- 
ings in which the overcrossing curve goes to the right of the inter- 
section are called positive. Those in which the overcrossing curve 
goes to the left of the intersection are called negative. The link- 
ing number is one-half the number of positive crossings minus the 
number of negative crossings. The linking number is symmetric, 
so that = 

The fact that this problem is #P-hard is a textbook exercise. See 
C. Moore and S. Mertens, The Nature of Computation, Ex. 13.11. 



